Formation and Stability of Synaptic Receptor Domains 
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Neurotransmitter receptor molecules, concentrated in postsynaptic domains along with scaffold 
and a number of other molecules, are key regulators of signal transmission across synapses. Employ- 
ing experiment and theory, we develop a quantitative description of synaptic receptor domains in 
terms of a reaction-diffusion model. We show that interactions between only receptor and scaffold 
molecules, together with the rapid diffusion of receptors on the cell membrane, are sufficient for the 
formation and stable characteristic size of synaptic receptor domains. Our work reconciles long-term 
stability of synaptic receptor domains with rapid turnover and diffusion of individual receptors. 



PACS numbers: 87.16.-b, 82.40.-g, 87.19.lp, 87.19.lw 

How the physiological stability necessary for memory 
storage can be achieved in the presence of rapid molec- 
ular turnover and diffusion is a central problem in neu- 
robiology [1 . Synapses, in particular, are believed to be 
the physiological seat of memory, and rely on the stabil- 
ity of postsynaptic domains containing neurotransmitter 
receptor molecules, as well as scaffold and a number of 
other molecules, over days, months, or even longer peri- 
ods of time ^ . Yet, recent experiments have demon- 
strated that individual receptor [4]-[6] and scaffold [7]- 
[9] molecules leave and enter postsynaptic domains on 
typical timescales as short as minutes. How can these 
seemingly contradictory observations — long-term stabil- 
ity and a well-defined characteristic size of postsynaptic 
domains on the one hand, rapid molecular turnover and 
diffusion on the other hand — be integrated in a unified 
understanding of postsynaptic domain formation and sta- 
bility? 

Classically, it has been assumed that interactions 
between presynaptic and postsynaptic neurons play a 
paramount role in the stability and in setting the char- 
acteristic size of synaptic receptor domains [10]. Over 
recent years, though, a number of studies [H HI [11] , car- 
ried out on a variety of chemical synapses, have indi- 
cated that molecular domains containing synaptic recep- 
tor molecules may form spontaneously even in the ab- 
sence of presynaptic neurons. However, a detailed molec- 
ular understanding of the mechanism governing the for- 
mation and stability of synaptic receptor domains has re- 
mained elusive. In this Letter, we first discuss a minimal 
experimental system which enables us to determine the 
molecular components essential for the self-organization 
of stable receptor domains of the characteristic size ob- 
served in neurons \9\ [121 (El- On this basis, we then 
formulate a mathematical model of the formation and 
stability of synaptic receptor domains which quantita- 
tively accounts for our experimental observations, and 



also makes further predictions pertaining to the stability 
and regulation of synaptic receptor domains. 

In our experiments we used single fibroblast cells, 
which are devoid of the molecular machinery commonly 
associated with postsynaptic domain formation [10 but 
allow for the rapid turnover and diffusion of receptors 
observed in neurons [H [6] , as well as for interaction of 
receptors with scaffold molecules. Fibroblast cells were 
transfected [14 with glycine receptors, one of the main 
receptor types at inhibitory synapses, and their associ- 
ated scaffolds, gephyrin molecules [9 . In our minimal 
system, the mere presence of both receptor and scaf- 
fold molecules led to the spontaneous emergence of sta- 
ble receptor-scaffold domains (RSDs) [see Fig. [TJa,b)]. 
These domains corresponded to a joint enhancement of 
the receptor and scaffold molecule densities, over a char- 
acteristic area of 0.2 to 0.3 /im^ [Fig. [ijc)]. Once the 
RSDs were formed, their mean area remained stable over 
a time scale of days, with little cell-to-cell variability in 
the mean area of RSDs but larger variability in the mean 
number of RSDs per cell [Fig. [ijc)]. If only receptors 
were transfected, in the absence of scaffold molecules, re- 
ceptor domains did not emerge, apart from possible oc- 
currences of transient microdomains [T2",T5^. If only scaf- 
fold molecules were transfected, in the absence of recep- 
tors, then these formed large intracellular blobs but no 
association with the cell membrane was detected [T2j [16] . 

The experiments carried out on our minimal system in- 
dicate that receptor-scaffold interactions, together with 
the diffusion properties of each molecular species at the 
membrane, are necessary and sufficient for RSD forma- 
tion and stability. In particular, the presence of a presy- 
naptic terminal is not essential for the occurrence of sta- 
ble RSDs. In agreement with previous studies f4l [6l [TI], 
our results point to a picture in which postsynaptic do- 
mains form in the absence of presynaptic stimulation, 
which subsequently intervenes in their maturation and 
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FIG. 1: (color online). Experimental results on the formation 
and stable characteristic size of RSDs composed of glycine 
receptors and gephyrin scaffolds [14 . (a) Example of a trans- 
fected COS-7 cell with domains on its membrane: Receptor 
(R, red) and scaffold (S, green) concentration and overlay (co- 
domains in yellow). A fraction of the apparent green- labeled 
scaffolds is endoplasmic. Scale bar, 5 /xm. (b) Examples of se- 
lected RSDs at higher resolution. For ease of visualization, the 
concentration maps of the two molecular species were slightly 
shifted with respect to one another in the color panel. Scale 
bars, 0.5 /xm. (c) Mean RSD (cluster) area and number of 
RSDs (clusters) per cell versus time. Error bars: standard er- 
rors; n > 10 cells from two independent experiments for each 
point. 



regulation. Both the characteristic size and the stability 
of the RSDs observed in our experiments are similar to 
those of synaptic receptor domains in neurons. Indeed, 
when scaffold molecules are transfected to young neu- 
rons devoid of synapses, domains of a comparable size 
arise jT2]. When they are transfected to mature neurons 
with synapses, the domain size remains unchanged [13]. 



Finally, the diffusion properties of receptors are similar 
in cells with transfected [9 and endogenous [17 scaffold 
molecules. Thus, we expect that receptors and scaffolds 
in neurons exhibit the necessary and sufficient proper- 
ties for RSD formation and stability, as they do in our 
experiments. 

We now turn to the mathematical description of our 
minimal experimental system. The concentration of re- 
ceptors is represented by the function r{x^y^t) and that 
of scaffolds by the function s{x,y,t)^ where the variables 
X and y denote coordinates along the cell membrane, and 
the variable t denotes time. The spatiotemporal evolu- 
tion of these fields is governed by the reaction-diffusion 
equations 



dr 

di 
ds_ 

di 



F{r, s) + i/^V [(1 - 5) Vr + rVs] , (1) 
G(r, s) + V [(1 - r)\/s + sVr] , (2) 



where F and G are simple (cubic) polynomials in r and 
s that describe the reactions in our system [14 , and 
and Us are the receptor and scaffold diffusion coefficients. 
The nonlinear corrections to the standard diffusion terms 
i/^V^r and z/gV^s in Eqs. ^ and ^ arise from the con- 
straint < r + 5 < 1 , where we have normalized r and s 
so that the maximum concentration of receptors and scaf- 
folds is equal to 1, and account for steric repulsion [4l[6] 
of receptors and scaffolds in the confined membrane envi- 
ronment of a living cell. Experimental studies jlHOllTS] of 
the diffusion properties of glycine receptors and gephyrin 
scaffolds, as well as of other types of synaptic receptors 
and scaffolds, yield Vr > J^s- 

The reaction and diffusion properties of receptors and 
scaffolds [IHSl [15] suggest that Eqs. ([T]) and ^ exhibit 
pattern formation via a Turing instability [18l |19] which 
emerges from the interplay between the two molecular 
species in our system: Receptors diffuse quickly and tend 
to repel nearby molecules, whereas scaffolds diffuse more 
slowly and tend to attract nearby molecules. In agree- 
ment with experiments, the formation of synaptic recep- 
tor domains via a Turing mechanism necessarily relies 
on the presence of both receptors and scaffolds. Expres- 
sions of the reaction terms F and G in Eqs. ([T]) and ([2| 
are obtained from the relevant chemical interactions, re- 
ported previously [4]-[9] , together with the general math- 
ematical constraints associated with Turing instabilities 
[19], a point we return to below. Reaction-diffusion mod- 
els akin to the one described here have, in recent years, 
been used to describe molecular localization during cell 
division [2QH22] . and are to be contrasted with models 
of domain formation which rely on phase separation and 
coarsening [23 . 

We simulated Eqs. ([T]) and ([2|, starting from random 
initial conditions, with units of space and time set by the 
diffusion coefficient of unbound receptor molecules and 
the rate of receptor endocytosis. Using typical values for 
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FIG. 2: (color online). Model results on the formation and 
stable characteristic size of RSDs 14 . (a) Irregular patterns 
of stable RSDs, with an area of approximately 0.2 to 0.3 /xm^ 
each, emerge on a timescale of hours. The distributions of 
receptors (upper panel) and scaffolds (middle panel) are over- 
layed (lower panel), and the domain shapes and patterns are 
stable, (b) Formation and (c) shape of RSDs at higher reso- 
lution. Scale bars, 0.5 /xm. 



these parameters taken from experiments [4]-[9l [TSJ, we 
found that irregular patterns of stable RSDs similar to 
experimental ones emerged over a timescale of hours [see 
Figs. [2](a) andjlja)]. Individual domains resulted from 
a coordinated increase of receptor and scaffold densities. 
Occasionally, we observed sets of closely-spaced RSDs in 
the outcomes of simulations, resulting from initial ran- 
dom fluctuations [see Figs. |2|a,b)]. At lower resolution, 
these appeared as larger than average, irregularly shaped 
domains reminiscent of similar instances obtained exper- 
imentally [see the lower panels in Fig. [ijb)]. Inclusion 
of molecular noise in our reaction-diffusion model is ex- 
pected to further distort the simulated patterns with- 
out changing their overall characteristics [24 . Moreover, 
we observed in our simulations that receptor aggrega- 



tion trails behind scaffold aggregation in time [Fig.|2|b)], 
which is in fact a general feature of our model and is also 
in agreement with experimental observations [25l [26] . 

The Turing mechanism implies that RSDs attain 
a steady state once a dynamical equilibrium between 
strong receptor diffusion and the attractive effect of scaf- 
folds is reached. In our simulations of Eqs. ([T]) and ^ 
we found that, for the diffusion coefficient of glycine re- 
ceptors reported in experiments [4]-[6l |9l [15] , the resulting 
characteristic size of RSDs was comparable to that ob- 
tained in experiments [see Figs. [2|b,c) and[TJc)]. This 
result necessitated only gross adjustment of the reaction 
rates, so that Eqs. ([T]) and ^ would exhibit a Turing 
instability but, apart from that, the characteristic size of 
domains was found to be largely insensitive to variations 
in reaction rates or, indeed, to the reaction kinetics con- 
sidered. The timescale for the formation of stable RSDs 
in our simulations [see Fig. |2|b)] was set by the rate of 
receptor endocytosis. For the measured range of values 
of the rate of receptor endocytosis [11 |6l |8] this timescale 
was also in broad agreement with experimental observa- 
tions [see Fig.Jllc)]. Further consistency checks between 
Eqs. ([T]) and ([2j and experimental results would pertain 
to the geometry of domain edges, the effect of changes in 
diffusion rates on domain size and stability, and the rel- 
ative extent of the enhanced concentrations of receptors 
and scaffolds in RSDs. In particular, our model predicts 
that once dynamical equilibrium has set in, the recep- 
tor profile is wider than the scaffold profile in any given 
domain [see Fig. |2|c)]. 

While, in accordance with experimental observations 
jlHOlfTS], we allowed [14 for a variety of interactions be- 
tween receptors and scaffolds when simulating Eqs. ([T]) 
and ([2|, we found that only a handful of chemical reac- 
tions were crucial for the formation and stability of RSDs 
via a Turing instability. To lowest order, these reactions 
correspond to R ^ Ri^ and R^ -\- S R -\- S for the 
receptors, and to 6* ^ ^'5 and Sb -\- 2S 3S for the 
scaffolds, respectively. In these expressions, the symbols 
R and S stand for receptors and scaffolds at the mem- 
brane, and Rb and Sb denote molecules in the bulk of the 
cell. In particular, the reaction -\-2S ^ 3S^ in which 
a scaffold molecule from the bulk is adsorbed onto the 
membrane by two other scaffold molecules into a trimer, 
is key to domain formation, whereas the simpler reaction 
Sb -\- S 2S alone is not sufficient. Indeed, gephyrin 
scaffold molecules are thought to form both dimers and 
trimers under the usual conditions in which neural do- 
mains are observed [4 . However, if trimerization is pre- 
vented, no domains (or only very small ones) appear [9]. 

The above results demonstrate how stable synaptic re- 
ceptor domains can emerge in the absence of presynaptic 
stimulation. In a synapse, however, presynaptic activ- 
ity regulates [27] the concentration of receptors in the 
postsynaptic domain. Our reaction-diffusion model sug- 
gests novel mechanisms for how such regulation may be 
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FIG. 3: (color online). Model results on the regulation of ma- 
ture RSDs J^. Left panel: temporal profiles of step stimula- 
tions. Inset: receptor concentration profiles at times of maxi- 
mum domain size. Right panel: time course of the in-domain 
receptor population size, R, following stimulation, normalized 
by the in-domain receptor population size in the absence of 
any stimulation, Rq. 

achieved. It has been observed [8l|28] that the diffusion of 
receptors on the postsynaptic membrane can be modified 
through binding of presynaptic neurotransmitters. Sim- 
ilarly, scaffold diffusion may be [29l [30] modulated by 
synaptic activity. This suggests that local modification 
of the diffusion properties of receptors or scaffolds may 
contribute to the regulation of postsynaptic domains. As 
a simple phenomenological perturbation to our model, 
we therefore implemented pre- and postsynaptic inter- 
actions through a local increase in the receptor diffu- 
sion rate which, within the framework of our reaction- 
diffusion model, has the same effect as a local decrease 
in the scaffold diffusion rate. 

As shown in Fig. [3) our model predicts that modu- 
lation of receptor or scaffold diffusion yields a transient 
increase of the in-domain receptor population following 
presynaptic activity. Clearly, this speculative, purely bio- 
physical mechanism may coexist with biochemical mech- 
anisms of postsynaptic plasticity. Applying a few seconds 
of stimulation at a time (Fig. [sj left panel), we found 
a correlation between the increase in in-domain recep- 
tor population and the duration of stimulation (Fig. [sj 
right panel). After an initial, short-lived suppression, the 
population increase lasted for a few tens of seconds — the 
timescale typically associated with short-term plasticity. 
In this window of time, RSDs were richer in receptors 
(Fig. [sj inset) and, hence, yielded a larger synaptic ef- 
ficacy. This phenomenon has a simple explanation in 
terms of the Turing instability exhibited by our model: 
Enhanced (diminished) receptor (scaffold) diffusion de- 
pletes the receptor (increases the scaffold) population in 
a transient manner which, because receptors are repulsive 
and scaffolds are attractive, in turn attracts even more 
receptors and scaffolds into RSDs. 

In summary, we have used a minimal experimental 
system devoid of synaptic machinery to show that neu- 
rotransmitter receptor domains of the stable character- 



istic size observed in neurons can emerge from noth- 
ing more than interactions between receptors and scaf- 
folds, together with the rapid diffusion of receptors on 
the cell membrane. A reaction-diffusion model quanti- 
tatively accounts for our experimental results, yielding 
spontaneous formation of stable receptor domains and 
their observed characteristic size, as well as new puta- 
tive mechanisms for the regulation of synaptic strength. 
Collectively, our results show how stable synaptic recep- 
tor domains may form even in the absence of presynaptic 
stimulation [H |6l [11] , and how rapid turnover and diffu- 
sion of receptors [4l|6], far from being a hindrance, may 
in fact be crucial [1] for ensuring overall stability and 
delicate control of synaptic receptor domains. 
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